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Universal quantum computers promise a dramatic speed-up over classical computers but a full- 
size realization remains challenging. However, intermediate quantum computational models have 
been proposed that are not universal, but can solve problems that are strongly believed to be 
classically hard. Aaronson and Arkhipov have shown that interference of single photons in random 
optical networks can solve the hard problem of sampling the bosonic output distribution which 
is directly connected to computing matrix permanents. Remarkably, this computation does not 
require measurement-based interactions or adaptive feed-forward techniques. Here we demonstrate 
this model of computation using high-quality laser-written integrated quantum networks that were 
designed to implement random unitary matrix transformations. We experimentally characterize the 
integrated devices using an in-situ reconstruction method and observe three-photon interference 
that leads to the boson-sampling output distribution. Our results set a benchmark for quantum 
computers, that hold the potential of outperforming conventional ones using only a few dozen 
photons and linear-optical elements. 



I. INTRODUCTION 

More than a decade ago, Knill, Laflamme and Milburn 
(KLM) showed in their seminal work [T] that scalable 
photonic quantum computing is possible using only lin- 
ear optical circuits, single-photon sources and detectors, 
and measurement-induced effective nonlinearities. The 
use of ancillary photons and adaptive feed-forward tech- 
niques [2] not only enables the heralding of successful gate 
operations IJ but also provides a basis for protocols in 
which probabilistic two-photon gates are teleported into 
a quantum circuit with high probability [S] . This has en- 
couraged researchers to pursue the realization of large- 
scale photonic quantum computers |6j. Even though 
there has been impressive theoretical progress [7HTU] , the 
required number of indistinguishable ancilla photons for 
a universal optical quantum computer |12j appears to be 
very challenging given current photonic quantum tech- 
nology. 

On the other hand, several interesting intermediate 
models of quantum computation have been recently pro- 
posed P^HTS] . Even though they do not enable universal 
quantum computation, these models still provide a dra- 
matic computational speed-up for particular tasks. Such 
quantum computers can be seen in analogy to quantum 
emulators that are designed for the simulation of one par- 
ticular physical system [TT]. In contrast to the proposed 
KLM scheme, these models need neither entangling gate 
operations, adaptive measurements, nor ancilla photons, 
and are thus technically more feasible. The intermediate 
quantum computation model proposed by Aaronson and 
Arkhipov jT6j seems to be extremely resource efficient as 
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it utilizes the unique advantages of the photons' mobility 
and bosonic nature to solve sampling problems that are 
beheved to be classically hard [171 IHI • 

Here we experimentally demonstrate the boson- 
sampling computation based on non-interacting and 
identical bosons. It is strongly believed that there is no 
known way of efficiently simulating the output of this 
computation by classical means [161. It is realized by 
using a quantum system composed of n non-interacting 
bosons that are processed through a network of m 
physical modes, where m > n. The bosonic nature of the 
photons themselves leads to non-classical interference 
when the photons propagate through a random network 
and produce a complex output probability distribution 
that is hard to sample on a classical computer. The 
simplest case of the non-classical interference effect was 
demonstrated for two photons in the seminal experiment 
of Hong-Ou-Mandel (HOM) [H] in 1987. When two 
single photons enter, one from each input mode, a beam 
splitter with 50% reflectivity (50/50 beam splitter), they 
will always exit the beam splitter together in one of the 
two output modes when the indistinguishability and 
temporal overlap is perfect (Figure la). This quantum 
physical phenomenon has not only become a standard 
technique for aligning quantum photonic devices, but 
also provides, together with the measurement process, 
the main underlying physical mechanism for non-trivial 
quantum gates in linear optical quantum computing |20) . 

This tendency of photons to "bunch" is explained 
by an effect in which two possible outcomes interfere 
destructively. Because a beam splitter represents a 
unitary transformation [21] . the probability that two 
photons entering a 50/50 beam splitter in input modes a 
and b respectively that both will be transmitted cancels 
with the probability that they will both be reflected. 
More generally, the probability of finding one photon in 
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output mode a' and the other in mode h' is given by the 
permanent of the beam sphtter matrix BS: 



P = |Pcr(BS)p = |Per 



T iR 
iR T 



|2 = |r2_^2|2 



(1) 



where T and R are beam sphtter transmission and re- 
flection coefficients, respectively. In the case of a 50/50 
beam splitter the permanent is obviously zero, therefore 
the photons bunch into one of the output modes. The 
same formula holds for the general case of n photons in- 
jected into n different modes of an m x m optical network 
with a underlying matrix U . The probability that one 
finds these photons in n different output modes is given 
by the permanent of the nxn sub-matrix U of the unitary 
U, P = |Per(t/)p (see Methods for details). For complex 
networks, such as randomly designed networks, comput- 
ing the permanent of the underlying unitary matrix on 
a classical computer is conjectured to scale exponentially 
in time [16j with respect to the size of the unitary ma- 
trix. The same holds for the output distribution. Thus, 
even for today's most powerful conventional computers, 
this puts an upper limit on the size of a unitary matrix 
for which the output distribution can be calculated. In 
particular, Rhode and Ralph [T7] estimated that for a 
random optical network of ti « 20 photons in m « 400 
modes, sampling the output distribution is already in- 
tractable for conventional computers. Obviously, these 
phyical requirements seem to be feasible in the not-too- 
distant future given the recent progress in photonic quan- 
tum technology |22) . This visualizes the importance of a 
benchmark computation with a lower number of photons 
in random optical networks where the classical verifica- 
tion is still possible. 

Even though bosons tend to bunch, also known as the 
"boson birthday paradox" [25] . for networks with a suf- 
ficiently large number of modes (m ^ \/^)^ the prob- 
ability of detecting n photons in n spatially separated 
modes as n-fold coincidences dominates. This remark- 
able feature reduces the technological requirements as no 
number-resolving detectors are needed for this interme- 
diate model of quantum computation and makes a full- 
fledged boson-sampling computation more feasible in the 
near future. Therefore, in our experiment we consider 
only those measurement outcomes where three photons 
are detected in three spatially separated modes as three- 
fold coincidences 



II. EXPERIMENT 

The boson-sampling computation is demonstrated for 
different randomly designed optical networks. The com- 
putation is initialized by the insertion of three indistin- 
guishable photons, one in each input mode, into an in- 
tegrated circuit with flve input and five output modes. 
After the propagation through the waveguide structures. 
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Figure 1. Non-classical interference, a) Basic two-photon 
interference. When two indistinguishable photons enter a 
50/50 beam splitter the two possible outputs, both photons 
being transmitted or reflected, are interfering destructively 
(highlighted in red). Therefore the probability of finding two 
photons in different output modes a and b' is zero. This is de- 
scribed mathematically as the permanent of the 50/50 beam 
splitter unitary, b) Example of an experimental three-photon 
interference. Three photons were injected into input modes 
3, 4, and 5 of one of our optical networks. We measured four- 
fold coincidence counts between a trigger detector and output 
modes 2, 4, and 5 of our device, while simultaneously delaying 
the input photons in mode 4, and 5. The Dip shows a clear 
signature of genuine three-photon interference. The error for 
the experimental data follows a Poissonian distribution of the 
measured counts. The shaded area represents the Gaussian 
fit including errors. 



the corresponding output distribution is recorded by de- 
tecting all possible three-photon coincidence measure- 
ments, where each photon is found in a different output 
mode. 

Each integrated circuit was fabricated with a direct 
laser writing technique [5S] and consists of five spa- 
tial modes that are coupled by eight beam splitters and 
eleven phase shifters (see Methods). To obtain different 
optical networks, two parameters were randomly varied 
during the fabrication process: the phases, by adjusting 
the relative path length between optical elements, and 
the beam splitter ratios, by tuning the evanescent cou- 
pling among the modes. The schematic drawing of these 
integrated optical networks is shown in Figure 2. The 
three-photon input state was generated via the process 
of spontaneous parametric down-conversion [26j . 

The photon source was aligned to emit the entangled 
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state |$+) = {\H}a\H}b + \V)a\V)b) /V^, where H and V 
denote horizontal and vertical polarization, respectively, 
and a and b correspond to the two spatial modes. By 
pumping with higher power (700mM^ cw-equivalent) also 
two photon pairs are emitted as a four-fold emission, 
while the even higher-order emission is kept low. To en- 
able a triggered three-photon emission, two photon pairs 
must be emitted simultaneously into spatial modes a and 
6, resulting in 

= {\HH)a\HH)b + \HV)a\HV)i, 

+ \VV)a\VV)b)/VS. 

These photons are guided to a state preparation stage 
utilizing two polarizing beam splitters (PBSl and PBS2) 
such that a successful detection event in the trigger mode 
a" heralds the generation of the states \H)a'\H)i,'\V)b" 
or IVF);,'. Postselection on a four-fold-coincidence, con- 
sisting of the trigger event and three detection events in 
the output modes of the circuit, ensures that three pho- 
tons entered the waveguide in separate spatial modes. A 
half-wave plate in mode b" introducing a 90° rotation is 
used to render the photons indistinguishable in polariza- 
tion. Using mating adapters, the three photons can be 
inserted in any combination of three input modes of the 
polarization-maintaining fiber-array that is butt-coupled 
to the integrated device. A schematic of the experimental 
setup is shown in Figure 3. 



III. RESULT 

The multi-photon interference on chip [27j is controlled 
by three adjustable delay lines to temporally overlap 
the photons. Scanning the temporal delays results in 
a three-photon HOM-dip that acts as a strong signature 
of the three photons' non-classical interference (Figure 1 
b). The underlying unitary operation of the integrated 
circuits was reconstructed by using an adaption of a re- 
cently proposed method [25]. For each optical network, 
the 19 independent parameters were fitted to the exper- 
imentally acquired 25 single-photon probabilities and 40 
two-photon visibilities (see Methods for details). Figure 
4 depicts the experimental data and theoretical predic- 
tions of the boson-sampling computation for two different 
randomly designed integrated circuits. For each experi- 
mental data point, four- fold coincidence events were de- 
tected for 20 hours. The output distributions are in good 
agreement with the theoretical values obtained from the 
reconstructed unitary matrices. 



IV. CONCLUSION 

Our experiment presents the first benchmark quan- 
tum computation on randomly designed optical networks 
showing a boson-sampling computation. This interme- 
diate model of quantum computation is of particular 





r 

{ 


127|jm[ 




127|jm 


3_ 


127|jm[ 


4_ 


127|jm|^ 


5_ 




10 cm 



b) 




Figure 2. The optical networks, a) Schematic drawing: 
The circuit consists of five input modes (1 to 5), five out- 
put modes (1' to 5'), eight directional couplers (771 to rjs) 
and eleven phase shifters (<^i to v'li)- Up to three single 
photons can be coherently launched into any combination of 
input modes. Each output mode is connected to a single- 
photon detector and coincidences are recorded with a home- 
built FPGA-logic. Neighboring modes are separated 127 fim 
and the chip exhibits a total length of 10cm. Three different 
optical networks written on the same chip were used in the 
experiment, b) Fluoresence image: In order to visualize the 
light evolution in the network, coherent laser light at a wave- 
length of 633nm is launched into input modes 2 to 4 of an 
optical network. Color centers are excited by the propagating 
beam and emit fiuorescent light at a wavelength of 650nm. 
The Fluoresence signal is directly proportional to the propa- 
gating light intensity. 



interest as the bosonic interference of photons in ran- 
dom networks is already hard to simulate on conven- 
tional computers. In contrast to universal models of 
photonic quantum computers that rely on ancilla pho- 
tons, measurement-induced interactions, and adaptive 
feed-forward techniques, the boson-sampling computa- 
tion requires only passive optical elements. This relaxes 
the physical requirements significantly such that a con- 
tinuous improvement of current multi-photon sources and 
detection efficiencies as well as reducing the losses in inte- 
grated circuits, might lead to quantum computations in 
regimes where classical verification is no longer possible 
in the near future. 
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Figure 3. Three photon probabilities and reconstructed unitary matrices. Experimentally measured (red on the left) 
and theoretical (blue on the right) three-photon probabilities for the ten possible output combinations. Photons were injected 
into input modes 3, 4, and 5 of device no. 1 a) and input modes 3, 4, and 5 of device no. 2 c). The errors for the probabilities 
follow a Poissonian distribution. Absolute value squared of the reconstructed matrix elements of b) optical network no.l and 
d) optical network no. 2. 



V. METHODS 



Boson sampling computation. This intermediate 
model of quantum computation involves a quantum sys- 
tem of n non-interacting bosons operating between m 
physical modes, where m > n. We define the computa- 
tion basis states |zi, 12, . . . , im) , where ii+i2 + - ■ - im = n. 
For example, the state |2, 1, 0, 1) is the state of two bosons 
in the first mode and one boson in the second and fourth 
mode. The total number of states in the computational 
basis D — (™^"~^) is exponentially large in to. Since 
there is no interaction, all the gates in the model are 
single-particle transformations defined by to x to com- 
plex unitary matrices. For the purposes of this work it is 
important to recall the definition of the permanent of a 



TO X TO matrix U 

m 

Per([/)= ^ H'^^-W' 

<jes„, 1=1 

where Sm is the set of all permutation of to elements. 
For example, the permanent of 2 x 2 unitary matrix U ~ 
icd) reads Per(J7) = ad + be. If the m x m matrix U 
exhibits sufficient complexity, e.g. is not built out of only 
and 1 elements, than the evaluation of its permanent 
is strongly believed to be hard on a classical computer. 
This means that the computation requires an exponential 
number of steps of computation in to. Consider an input 
state |/) = 12, . . . , im) of n bosons in to modes that 
is transformed via some unitary m x m matrix U. The 
probability of finding the state \0) = J2, ■ • • ,jm) at 
the output is given by 

Pjo = \{0\U(E}U---(E)U\I)\'^. 
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Figure 4. Experimental setup. The experimental setup can be divided into three main parts. The first part is a down- 
conversion source pumped with 150/s-pulses at a wavelength of 394. 5nm and cw-equivalent power of 700mW. Four photons 
are created using higher-order emissions from a 2mm BBO crystal cut for a non-coUinear typell process. Filtering with 
^FWHM ~ 3nm filters guarantees spectral indistinguishability. The photons are coupled to single-mode fibers and guided to 
the state preparation stage. Here, four-fold events are split up via two polarizing beam splitters (PBS) and a half-wave plate 
(HWP) in mode b" ensures that all photons entering the integrated device exhibit the same polarization. The input photons 
are coupled into polarization maintaining fibers (PM-fibers) , which can be mated to any of the five fibers of the input-fiber 
array. Temporal overlap in the circuit is achieved with three delay lines. After passing through the waveguide, the photons are 
coupled to graded-index multimode fibers and sent to single-photon avalanche photo diode detectors (Dl to D5). 



Careful analysis [TO shows that this probability can be 
expressed trough the matrix permanents: 



Pi,o - 



|Per(C//,o)P 



ii\i2\ 



where the n x n matrix J7/,o is defined in the following 
manner. First we define the mxn matrix Uj by taking 
copies of the fc'^ column of U, for each k = 1 . . .m. Then 
we form the n x n matrix ?7/.o by taking jk copies of 
the fc"^ row of Uj , for each k — 1 . . .m. As an example, 
consider a 3 x 3 matrix 




and the input \I) = 1 1,1,0) and output state \0) 
|0, 1, 1). Than the matrix Ui is 




and finally C//_o reads 



Uin = 



d e 
9 h 



Therefore, the probability of finding the state \0) = 
|0,1,1) at the output, that is to find one boson in the 
mode 2 and one in mode 3 is given by 



P; 



|Per(C/^ 



i,o) 



1!1!0!0!1!1! 



\dh — eg\ 



From the previous example it is clear that for cases where 
the input and output modes are occupied by at most one 
boson the matrix Ui^o turns out to be the n x n subma- 
trix of U. 

Experimental setup. An 80MHz Ti:Sapphire oscilla- 
tor (Chameleon, Coherent Inc.) emitting 150/s pulses 
at 789nm (2.5W cw-equivalent power) is upconverted to 
700mW cw-equivalent power at 349. 5nm via a LiBaOs 
crystal (LBO) (HarmoniXX, A.P.E GmbH). The beam 
is focused on a 2mm /3-BaB204 (BBO) crystal cut for 
degenerate non-coUinear type-II spontaneous parametric 
downconversion j26j . To achieve spectral indistinguisha- 
bility, the downconverted photons are filtered by interfer- 
ence filters {Xfwhm — 3rtm) and collected with single 
mode fibers. The cw-equivalent pump power of TOOmW 
allows for emission of four-fold states at high count rates. 
Noise contribution from higher-order terms was measured 
to be lower than four percent. The double-pair emission 
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of the source generates the foUowing state 



\^)a,b = {\HH)a\HH), + \VV)a\VV)t 
+ \HV)a\HV)b)/V3. 

Two polarizing beam spUtters (PBSl and PBS2) dis- 
tribute the photons into four modes such that a four- 
fold coincidence postselection guarantees that three pho- 
tons entered the waveguide in different modes. Half-wave 
plates ensure indistinguishability in polarization. Three 
delay lines are used to temporally overlap the photons 
for non-classical interference. The photons reach the in- 
tegrated circuit via a polarization maintaining v-groove 
fiber array that is butt-coupled to the waveguide. In- 
dex matching gel is applied to reduce reflection-losses. A 
graded-index multimode fiber array is butt-coupled to 
the output of the waveguide and connected to single- 
photon avalanche photo detectors (SASPDs). Optimal 
coupling between the integrated device and the fiber ar- 
rays is achieved with two six-axis alignment stages. 



ACKNOWLEDGMENTS 

The authors are grateful to S. Aaronson, C. Brukner 
and M. Ringbauer for discussions. We acknowledge sup- 
port from the European Commission, Q-ESSENCE (No 
248095), QuILMI (No 295293) and the ERA-Net CHIST- 
ERA project QUASAR, the German Ministry of Educa- 
tion and Research (Center for Innovation Competence 
program, grant 03Z1IIN31), the John Templeton Foun- 
dation, the Austrian Nano-initiative NAP PlatonAus- 
trian Science Fund (FWF): (SFB-FOCUS) and (Y585- 
N20) and the doctoral programme CoQuS, and the Air 
Force Ofhce of Scientific Research, Air Force Material 
Command, United States Air Force, under grant number 
FA8655-11-1-3004. 



Chip fabrication. The waveguides were written inside 
high purity fused silica (Corning 7980ArF grade) using 
a RegA 9000 seeded by a Mira Ti:Sapphire femtosecond 
laser. The pulse duration was 150/s at 800nm with a 
repetition rate of lOOkHz and a pulse energy of 200?! J. 
The pulses were focused 370^m under the sample surface 
using a NA = 0.6 objective while the probe was trans- 
lated with a constant speed of 6cm /min by high-precision 
bearing stages (ALS130, Aerotech Inc.). The mode field 
diameter of the guided modes is 22^m x 22fj,m at 789nm. 
Propagation loss was measured at Q.3dBcm~^ and bire- 
fringence was measured in the order of i? = lO"*". Mea- 
sured coupling with the input fibers (850nTO PM-fibers, 
OZ Optics) was —6.7dB while losses at the output facet 
(graded-index multimode with 50/zm core diameter) were 
negligible. 

Unitary matrix reconstruction. To verify the exper- 
imentally obtained output distributions with theoretical 
data, the random 5x5 unitary matrices had to be recon- 
struced. Extracting the actual unitary transformation 
experimentally |28H30j is important as small imperfec- 
tions in the fabrication process can lead to significant 
changes in the transformation matrix. Since the position 
of beam splitters and phase shifters were known for our 
integrated devices and the output loss was experimen- 
tally extracted, the fitting of the 19 unknown variables 
(eight beam splitter ratios and eleven phase shifts) could 
be done via single-photon data and two-photon HOM- 
dips [28] . All errors are conservatively computed by as- 
suming Poissonian error propagation. 
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